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Stable  Local  Nodal  Bases  for  C1 
Bivariate  Polynomial  Splines 


Oleg  Davydov  and  Larry  L.  Schumaker 


Abstract.  We  give  a stable  construction  of  local  nodal  bases  for  spaces 
of  C 1 bivariate  polynomial  splines  of  degree  d > 5 defined  on  arbitrary  tri- 
angulations. The  bases  given  here  differ  from  recently  constructed  locally 
linearly  independent  bases,  and  in  fact  we  show  that  stability  and  local 
linear  independence  cannot  be  achieved  simultaneously. 


§1.  Introduction 

Given  a regular  triangulation  A,  let 

Sj( A)  :=  {s  £ Cr(ft)  : s|t  € Vd  for  all  triangles  T £ A}, 

where  Vd  is  the  space  of  polynomials  of  degree  d,  and  ft  is  the  union  of  the 
triangles  in  A.  In  this  paper  we  focus  on  the  case  r = 1 and  d > 5.  The  main 
result  of  the  paper  is  a construction  of  a basis  B {I?j}"=1  for  <Sj(A)  with 
the  following  properties: 

PI)  The  basis  B is  local  in  the  sense  that  for  each  1 < i < n,  the  support  of 
Bi  is  contained  in  star(uj)  (see  the  end  of  this  section)  for  some  vertex  Uj, 
P2)  The  set  B is  stable  in  the  sense  that  there  exist  constants  K\  and  K2 
dependent  only  on  the  smallest  angle  9&  in  A such  that 

n 

tflMoo  < llI>.EM|oo  < A~2  ||c||oo  (1) 

t=l 

for  all  choices  of  the  coefficient  vector  c = (ci, . . . , cn). 

Bases  for  <S](A)  satisfying  property  PI  were  constructed  in  [14]  using  nodal 
techniques,  but  they  fail  to  satisfy  property  P2  for  triangulations  with  near 
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singular  vertices  or  near  degenerate  edges,  even  if  the  smallest  angle  in  the 
triangulation  is  controlled. 

For  convenience,  we  recall  the  definitions  of  some  of  the  terminology  used 
above.  Suppose  v is  a vertex  of  a triangulation  which  is  connected  to  Vi , V2 , V3 
in  counter-clockwise  order.  Then  the  edge  e :=  (v,  112)  is  said  to  be  near- 
degenerate at  v (degenerate  at  v)  provided  that  the  edges  (v,vi)  and  {v,V3) 
are  near-collinear  (collinear).  The  vertex  v is  called  near-singular  (singular) 
if  there  are  exactly  four  near-degenerate  (degenerate)  edges  attached  to  it. 
Given  a vertex  v of  A,  star(u)  = star^u)  is  the  set  of  triangles  sharing  v, 
and  starf(i>)  is  defined  recursively  as  the  union  of  the  stars  of  the  vertices  of 
star£-1(v). 


§2.  Nodal  Determining  Sets  and  Nodal  Bases 

Suppose  s is  a spline  in  Sj( A),  and  that  v is  a point  in  fi.  In  this  paper 
we  are  interested  in  certain  linear  functionals  defined  on  <Sj  ( A)  in  terms  of 
values  and  derivatives  of  s at  points  v in  ft.  Such  functionals  are  called  nodal 
functionals.  There  are  three  types  of  nodal  functionals  of  interest  here: 

1)  the  value  s(v), 

2)  the  directional  derivative  D™s(v),  where  w is  a given  vector  and  m is  a 
positive  integer, 

3)  the  mixed  derivative  D™1  D™£s(v)  at  a vertex  v of  A,  where  wi  and  w 2 
are  two  noncollinear  vectors  which  point  into  a common  triangle  T of  A. 

Definition  1.  A collection  M :=  (Aj}"=1  of  nodal  functionals  is  called  a 
minimal  nodal  determining  set  for  <S](A)  provided  they  form  a basis  for  the 
dual  space  (<S^(A))*.  If  M is  such  a set,  then  there  exist  unique  splines 
B in  S](A)  such  that 

\{Bj  — 6ij , i,  j — 1, . . . , 71.  (2) 

We  call  B a nodal  basis  for  S^( A). 

In  this  paper  we  will  concentrate  on  nodal  functionals  which  involve 
derivatives  De  along  edges  e :=  (yi , V2)  of  the  triangulation  A,  or  perpen- 
dicular to  such  edges.  Denoting  the  Cartesian  coordinates  of  a point  v by 
(vx,vy),  we  see  that  the  derivative  along  the  edge  e is  given  by 

n , (nf  - vf)Dxs(v)  + ( vv2  - vl)Dys{v) 

V(vx2-vxi)2  + (vy2-vl)2  ' 

while  the  derivative  perpendicular  to  the  edge  e is  given  by 

j~\  (^2  - vVi)Dxs{v)  - (vf  - ^)ZV(v) 

UpLSyv ) :=  , ■ = . 

v/K  - vfy  + (vy  - 

Note  that 

D(vi,V2)S(V)  — ~£*(v2,Vi)S(V)t  ^{v1,V2)-LS(V^  = ~^{V2,Vi)±S(V)' 
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§3.  Smoothness  Conditions  Between  Polynomial  Pieces 

It  is  well-known  how  to  describe  smoothness  between  polynomials  defined  on 
adjoining  triangles  in  terms  of  the  Bernstein-Bezier  coefficients  of  the  two 
polynomials.  Here  we  need  similar  conditions  in  terms  of  nodal  information. 
Suppose  T = (vi,V2,vz)  and  T = (vi,v2,V3)  are  two  adjacent  triangles  with  a 
common  edge  e = (vi,v2).  We  set  8\  = lv$v  iv2,  O2  = Lv^v^vi,  6\  = Lv-iv\V2, 
62  — /.V3V2V1.  Suppose 


vi  < v{'°  < ■■■  Vj’°5  < v2 

(3) 

vi  < vl’1  < ■■■  < v2 

are  given  points  lying  in  the  interior  of  the  edge  e. 

Lemma  2.  Letj),  p be  polynomials  of  degree  d > 5 defined  on  adjoining 
triangles  T and  T as  above.  Then  p and  p join  together  with  smoothness  C1 
across  the  edge  e :=  (vi,v2)  if  and  only  if  the  difference  g = p — p satisfies 


g(vi)  = Deg(vi)  = De±g(vi ) = D*g(vi)  = 0, , i = 1, 2, 

g(vi'°)  = °>  * = 1, . . . ,d  — 5, 
De-'-di'Vi’1)  = * = 1, . . . ,d  — 4, 

and 


>p(vl)  — &lD(VltV2)D(VltV3)p{vi)  + <JiD{VltV2)D{Vlyj3)p(vi), 
°2D(v2,vi)P(v2)  = ^2D{V2tVl)D^V3)p{v2)  + a2D^Vx)D^i3)p{v2), 


(4) 

(5) 

(6) 


where  o,  :=  sin  0l,  d,  :=  sin  6,,  d,-  :=  sin(d,  + 9i),  i = 1,2. 

Proof:  We  follow  the  method  of  proof  of  the  main  result  in  [14].  Concerning 
necessity,  we  first  observe  that  if  p and  p join  with  C1  continuity  across  e, 
then 

g(v)  = Dwg(v)  = 0,  for  all  v € e,  (7) 

where  w is  any  unit  vector  noncollinear  with  the  edge  e.  This  implies  (5) 
and  the  conditions  on  g,  Deg  and  De±g  in  (4).  The  conditions  on  the  second 
derivatives  are  easily  obtained  by  differentiating  the  identities  (7)  along  the 
edge  e and  using  the  fact  that 

diD( 

Vl,T>2  )P{v  1)  = di  D(  Vl,V3  )P(ui)  + ffifi(t,i,i3)p(r1), 

d2T>( 

V2,Vi  )P{v  2)  = d2T>(  V2,V3  )p{v2)  + a2D{v2,v3)P{v2). 

To  prove  sufficiency,  suppose  that  p and  p satisfy  (4)-(6).  Then  the 
univariate  polynomial  g\e  is  of  degree  at  most  d and  satisfies  d + 1 homoge- 
neous Hermite  interpolation  conditions  on  e.  Therefore  g(v)  = 0 for  v £ e. 
This  shows  that  p and  p join  continuously.  We  now  consider  the  cross- 
derivative q :=  Dexg\e  which  is  a univariate  polynomial  of  degree  at  most 
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d—  1.  By  (4)— (5) , q has  d— 2 zeros  . . . ,ve^\,V2  on  e.  Moreover,  by  (6), 

Deq(v i)  = Deq{v 2)  = 0,  as  is  easy  to  check  by  expressing  DeDe±p(vi)  as  a lin- 
ear combination  of  D%p(v  1)  and  DeD^VlV3^p(vi)  and  expressing  DeDe±p(vi) 
in  terms  of  Dlp(v  1)  and  DeD^Vltg3yp{v  1),  and  similarly  for  V2 ■ Therefore, 
q = 0,  and  we  have  shown  that  p and  p join  with  C'1-smoothness.  □ 

For  a different  set  of  nodal  smoothness  conditions,  see  [5]. 


§4.  Construction  of  a Stable  Local  Nodal  Basis  for  A) 

In  this  section  we  begin  by  defining  a spanning  set  A/a  of  nodal  functionals 
for  (<5>](A))*.  Then  we  choose  an  appropriate  linearly  independent  subset 
M which  forms  a basis  for  (S]( A))*.  This  will  involve  analysing  the  linear 
dependencies  between  elements  of  A/a  (he-,  the  smoothness  conditions).  The 
corresponding  nodal  basis  determined  by  the  duality  conditions  (2)  will  be  the 
desired  stable  local  basis  for  <Sj( A).  Given  a triangle  T :=  (rq,  V2,  ^3),  let 
iv  1 + jv  2 + kv3 

vijk  j 1 


i + j + k = d. 


Given  an  edge  e of  A,  let  v®’°  and  v be  the  points  defined  in  (3).  We  define 


nT  ._  r \T  „ _ 


s(vTjk)  '•  i + j + k = d , 2 <i,j,k}, 


£{e)  ~ {\fs  = s{vf)  : i = 1, ...  ,d  — 5} 

U{A*’1s  = \e\DeA.s{Vi'1)  : i = 1, . . . , d - 4}, 
where  |e|  denotes  the  length  of  e. 

Given  a vertex  v in  A,  suppose  the  vertices  connected  to  v are  iq , . . . , vn  in 
counterclockwise  order  (with  rq  a boundary  vertex  if  v lies  on  the  boundary), 
and  let  Tl'l  = (v,Vi,Vi+i),  e*  = ( n,rq ),  9{  = Zejei+i,  where  if  v is  an  interior 
vertex,  we  identify  vg+n  = vg,  eg+n  = eg.  Denote  by  |star  (v)\  the  diameter  of 
star(n).  Let 

T>i(v)  :=  {A?s  = |star(n)|,+-'DJ.jC>j)s(t;)  : 0 < i + j < 1} 

■£•2  W :=  {A?ps  = sing.Sin8i_,  |star  {v)\ 2D2Cis{v)  : i = 1, . . . , n} 

U{A-,ms  = jjAg-|star (n)|2Dei Dei+1s(v)  : i = l,...,n] 
if  v is  an  interior  vertex,  and 


Mv)  ■■=  (A^p 


rw47rc-|star(v)|2D2.s(v)  : i = 2,. ..  ,n  - 1} 


u {Ar,Ps  = |star(v)|2D2.s(?;)  : i = 

U{A”ms  = jj^jrlstar (v)\2 De.Dei+1s(v)  : i = 1} 

if  v is  a boundary  vertex.  Let 

A/a  :=  1J  CT{J  U £(e)u  U [PiOOUttaH]. 

Te  A eCA  veA 

The  Markov  inequality  implies  that  for  all  s € iSj(A)  and  all  A € A/a, 

|As|  < K\\S\U 

for  some  constant  depending  only  on  d and  the  smallest  angle  6 a in  A. 


(8) 
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Lemma  3.  The  set  A/a  is  a spanning  set  for  (5] (A))*.  Moreover,  the  only 
linear  dependencies  between  elements  of  A/a  are  given  by 

^i,m  + ~ sin(&i  + 0i-l)^i,p  (9) 

for  every  vertex  v and  every  interior  edge  e,  attached  to  v. 

Proof:  Let  s £ <S](A).  If  As  = 0 for  all  A € A/a,  then  on  each  triangle  T € A 
there  are  exactly  (d^2)  homogeneous  Hermite  interpolation  conditions  on  s, 
and  it  is  easy  to  see  that  they  force  s to  be  zero.  It  follows  that  A/a  is  a 
spanning  set  for  (<Sj( A))*.  The  second  statement  follows  immediately  from 
Lemma  2.  □ 

Algorithm  4.  (Construction  of  a stable  local  nodal  basis  for  <Sj(A).j  Let 
{Bi}"=1  be  the  set  of  splines  determined  by  the  duality  conditions  (2)  corre- 
sponding to  the  following  set  M :=  {A;}f=1  of  nodal  functionals: 

1)  For  each  triangle  T,  choose  the  (d~4)  nodal  functionals  Afjk  in  CT . 

2)  For  each  edge  e = (vi,  V2),  choose  the  2d  — 9 nodal  functionals  A®’0  and 

K’1  in  £(e)- 

3)  For  each  vertex  v,  choose  the  three  nodal  functionals  Ay  in  T>i(v). 

4)  For  each  vertex  v,  choose  the  following  nodal  functionals  in  72.2(f): 

a)  one  of  the  functionals  A”m  corresponding  to  the  first  mixed  derivative 
at  v,  and 

b)  all  functionals  Avi  p corresponding  to  the  pure  second  derivatives  at 
v,  with  one  exception:  if  v is  a nonsingular  interior  vertex,  the  func- 
tional A?  is  omitted,  where  io  is  chosen  such  that 

|sin(0io  + 0jo-i)|  > |sin(0j  + for  all  i = 1, . . . ,n.  (10) 

Theorem  5.  The  set  M.  of  Algorithm  4 is  a minimal  nodal  determining  set 
for  <S]( A),  and  the  nodal  basis  {Bi, . . . , B^}  for  <Sj(A)  defined  in  (2)  is  local 
and  stable,  i.e.,  it  satisfies  both  conditions  PI  and  P2. 

Proof:  The  fact  that  Ad  is  a basis  for  (Sj(A))*  follows  easily  from  Lemma  3. 
To  construct  a typical  basis  spline  Bj,  we  set  A {Bj  = <5y  for  all  i = 1, . . . ,n. 
Then  the  remaining  nodal  values  A Bj , A £ A/a  \ A4  are  computed  from  the 
smoothness  conditions  (9).  It  is  easy  to  see  that  the  support  of  the  resulting 
spline  is  at  most  the  star  of  a vertex.  This  shows  that  PI  is  satisfied. 

It  remains  to  show  that  the  Bj  form  a stable  basis.  This  follows  from  (8) 
by  a standard  argument  [12],  provided  we  can  show  that 

Halloo  < A,  l<j<n,  (11) 

where  A is  a constant  depending  only  on  d and  the  smallest  angle  9&  in  A. 
This  clearly  holds  if 


\ABj\<K,  for  all  A £ A/a 
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for  a similar  constant  K.  By  construction,  |ABj|  < 1 for  all  A € Ad.  Since 
A/a  \ Ad  C U«  ^{v),  let  us  take  an  arbitrary  vertex  v of  A and  notice  that 
if  A € 7Z2(v),  then  A Bj  can  be  nonzero  only  if  the  corresponding  Xj  lies  in 
77.2  (i>).  Therefore,  it  will  be  sufficient  to  show  that  | A B j | < K for  all  j such 
that  A j e 72-2  (v)  and  all  X elZ2{v)\M.  We  distinguish  four  cases. 

Case  1:  ( v is  a boundary  vertex.)  In  this  case,  ^(v)  \ Ad  = {A"m  : 
i = 1,  j / i i},  where  A?  m is  the  functional  included  in  Ad  in  step 

4a)  of  Algorithm  4.  Without  loss  of  generality  we  assume  that  ij  = 1.  For 
any  s g S\( A),  (9)  implies 

^2  ,ms  = ~K,ms  +<72^2,ps 

A3 ,ms  = A” )JnS  — 02A 2iPS  + CT3A3  pS 


A” 


n— 1 ,ma 


(-i)"ay>to«  + 

i=2 


where  we  set 

<7;  :=  sin(0i  + 0i_i). 

If  we  take  s to  be  the  basis  spline  Bj  corresponding  to  a A j € 722  (v),  then 
all  but  one  of  the  values  on  the  right-hand  side  of  the  expression  for  X\mBj 
vanishes,  and  thus 

lA^B.I  < |AjBj|  = 1,  * = 2, . . . ,n  — 1, 


which  proves  the  assertion. 

Case  2:  ( v is  an  interior  vertex  with  n / 4.)  In  this  case,  72-2 {v)  \ Ad  = 
{AVm  : i = i ± ij)  U {AV  p}.  For  AVms,  i = 1, . . . ,n,  i f h,  the 

same  calculation  as  in  Case  1 applies:  we  start  from  A”  ms  and  calculate  A)'ms 
consecutively  counterclockwise  until  AJ’o_1  ms,  and  then  also  clockwise  until 
A i0,ms-  For  A-0,ps,  we  have  by  (9), 


A’o^  = ffr„1(A?„,ra*  + A?0_1,m5). 


(12) 


Therefore,  our  claim  will  be  established  if  we  show  that 


k.o1!  = lsin  1i0io  +^0-1)1  < K3  if  n/4,  (13) 

where  K3  is  a constant  dependent  only  on  This  is  obvious  for  n = 3. 
Assuming  n > 5,  we  have  \6\  + 82  + #3  + 64  — 2ir\  > 6 a-  Hence, 


\6l0  + 0io-i  - 7r|  > max{|0j  +e2-  7r[,  1 63  +64-  tt|}  > 0A/2, 


and  (13)  follows. 
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Case  3:  (v  is  a singular  vertex.)  In  this  case,  Tl2{v)\M  = {A}m  : i = 2, 3,4} 
(where  we  assume  for  simplicity  that  *i  = 1).  Since  ay  = • ■ • = 04  = 0 for  a 
singular  vertex,  (9)  now  reduces  to 


K,ms  + — 0,  i — 1,2,3,  4. 

Therefore, 

Kms  = (-l)i+Xms-  * = 2,3,4, 
and  the  assertion  follows. 

Case  4:  ( v is  a nonsingular  interior  vertex  with  n = 4.)  We  proceed  as  in 
Case  2,  but  calculate  A}0  ps  differently.  At  first  glance  it  may  seem  that  (10) 
does  not  guarantee  stability  since  |crio  | may  be  arbitrary  small  (in  the  case 
of  near-singularity),  while  A"o  ps  is  to  be  computed  from  (12).  However,  the 
complete  system  of  equations  (9)  for  'R.^v)  is 


<7iAVps  = AVms  + AV_lims,  * = 1,2, 3, 4. 

Taking  the  sum  with  alternating  signs,  we  get 

D-i^a^o, 

i=l 

and  hence 

K,p°\  < E s 1 

i*io  l<Tto1 

for  every  s = Bj , with  A j € 7^2 (u).  This  completes  the  proof  of  (11),  and  the 
theorem  has  been  established.  □ 


§5.  Stability  vs.  LLI 

We  recall  (cf.  [2, 4, 6, 8, 9])  that  a set  B of  basis  splines  in  5}  (A)  is  called  locally 
linearly  independent  (LLI)  provided  that  for  every  T G A,  the  splines  {Bi  : i € 
Ex}  are  linearly  independent  on  T,  where 

Ex  :=  {*  : TCsuppHj}.  (14) 

A star-supported  LLI  nodal  basis  was  constructed  for  <Sj(A)  in  [4],  We  now 
establish  the  following  surprising  result. 

Theorem  6.  For  d > 5,  it  is  impossible  to  construct  a basis  for  <S}( A)  which 
satisfies  both  conditions  P2  and  (14)  simultaneously. 

Proof:  Suppose  [B\, . . . , Bn}  is  a locally  linearly  independent  basis  for  «S}(  A) 
on  a triangulation  A which  contains  an  interior  near-singular  vertex.  Suppose 
v is  connected  to  tq , . . . , V4  in  counterclockwise  order,  and  let  e,  be  the  edge 
(*y Vi),  Ti  the  triangle  (vi,vi+i,  v),  and  9j  the  angle  between  e,  and  e,+i- 
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Suppose  that  none  of  e,  is  degenerate  at  v.  For  each  1 < j < 4,  let  Sj  be  the 
unique  spline  in  Sj  (A)  such  that 


X i,msj  ~~  ^ij  ) b j — 1,  2,  3,  4, 

A Sj  = 0,  for  all  A £ A/a  \ 77-2  (v) 

Clearly, 

supp  sj  = Tj_i  U Tj  U Tj+ 1, 
and  we  can  write  (see  [2,9]) 

sj  = 


ieij 


where  Ij  :=  {z  : supp  fi;  C supps^}.  We  now  consider  the  spline 
s = -si  + s2  - s3  + s4  = - ^2  c^Bi  + ci21-8'  _ ^2  c^Bi  + ^2  c^Bi 

tG/i  i£h  *€/3  iGA 

J2  OfBi. 


ieiiuhuhuh 


Using  the  smoothness  conditions  (9),  it  is  easy  to  see  that 

A*pi  = 0,  AVmi  = (-l)i,  * = 1,2, 3, 4, 

As  = 0,  for  all  A € A/a  \ ^(v). 

and  thus  ||s||oo  < K 4,  where  K 4 depends  only  on  d.  If  the  basis  {f?i, . . . ,Bn} 
satisfies  P2,  we  get 

|Hloo<0|s|loo<tf4/tfl- 

Moreover,  since  A 2,mBi  7^  0 only  if  7i  U T2  U T3  C supp  B{,  we  have 


1 = ^2 ,m«  = aiA2 ,mBi  < #h\\a\\oo  max  \X%mBi\, 
ieh 

with  J2  :=  {z  : suppB,  = T\  U T2  UT3}.  Clearly,  #f2  < 3 (d^2) , and  hence 
there  exists  i 0 € I2  such  that 


\XlmBl0\  >K 5 > 0, 

where  FTs  depends  only  on  0a  ■ However,  X\  jnBUl  = 0,  so  that  by  (9)  we  have 


|A5. 


?Bio1  |sin(01  + 02)! 


IA2,, 


>Bi0 1 > 


K5 

sin(0i  + 02)|  ’ 


which  is  unbounded  as  0i  + 02  — > 77  ■ In  view  of  the  Markov  inequality,  it 
follows  that  HSJloo  is  unbounded.  But  then  the  basis  {B 1, . . . ,Bn}  cannot 
be  stable.  □ 


Stable  Local  Nodal  Bases 


179 


§6.  Remarks 

Remark  1.  Stable  local  bases  are  important  for  both  theoretical  and  practical 
purposes.  For  example,  it  can  be  shown  (see  [12])  that  if  a spline  space  has 
such  a basis,  then  it  has  full  approximation  power.  Applications  where  stable 
bases  are  useful  include  data  fitting  and  the  numerical  solution  of  boundary- 
value  problems. 

Remark  2.  For  d > 5,  stable  local  bases  for  certain  superspline  subspaces  of 
<S^(A),  can  be  constructed  using  classical  finite  elements,  see  [15].  However, 
it  is  also  important  to  have  such  bases  for  the  full  spaces  <S^(A),  since  in 
contrast  to  supersplines,  they  are  nested,  i.e.,  5^(Ai)  C <Sj(A2)  whenever 
A2  is  a refinement  of  Ai.  This  is  important  for  multiresolution  applications, 
see  [3,13]. 

Remark  3.  Algorithm  4 is  a modification  of  the  algorithm  used  in  [14]  to 
construct  a star-supported  basis  for  (A).  The  only  change  is  in  the  choice 
of  nodal  functionals  in  step  4b)  where  iq  was  taken  to  be  any  index  such  that 
eio  is  nondegenerate  at  v.  To  get  stability,  we  have  to  choose  »o  more  carefully. 
The  basis  constructed  in  Algorithm  4 is  not  locally  linearly  independent.  To 
get  an  LLI  basis,  step  4)  has  to  be  modified  in  a different  way,  see  [4]. 

Remark  4.  Star-supported  bases  were  constructed  for  general  spline  spaces 
S%( A)  for  d > 4r+l  in  [1],  and  for  d > 3r-f  2 in  [10,11].  The  constructions  were 
based  on  Bernstein-Bezier  techniques,  and  are  not  stable  for  triangulations 
that  contain  near-degenerate  edges  and/or  near-singular  vertices. 

Remark  5.  In  [7]  we  use  Bernstein-Bezier  techniques  to  construct  stable  local 
bases  for  general  spline  spaces  <SJ( A)  and  their  superspline  subspaces  for  all 
d > 3r  + 2.  In  a related  work  [6] , we  also  used  Bernstein-Bezier  techniques 
to  construct  locally  linearly  independent  bases  for  the  same  range  of  spline 
spaces  and  superspline  spaces.  For  more  on  LLI  spaces,  including  applications 
to  almost  interpolation,  see  [2, 4, 6, 8, 9]. 

Remark  6.  Following  the  arguments  in  [7],  it  is  easy  to  show  that  a natural 
renorming  of  our  stable  bases  is  Lp-stable  for  all  p £ [l,oo]. 
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